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We demonstrate that the dynamics of an open quantum system can be calculated efficiently and 
with predefined error, provided a basis exists in which the system-environment interactions are local 
and hence obey the Lieb-Robinson bound. We show that this assumption can generally be made. 
Defining a dynamical renormalization group transformation, we obtain an effective Hamiltonian for 
the full system plus environment that comprises only those environmental degrees of freedom that 
are within the effective light cone of the system. The reduced system dynamics can therefore be 
simulated with a computational effort that scales at most polynomially in the interaction time and 
the size of the effective light cone. Our results hold for generic environments consisting of either 
discrete or continuous degrees of freedom. 
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Introduction The interaction with its environment 
causes a quantum system to loose energy and phase - 
this is termed decoherence [T]. Since a quantum system 
can never completely be isolated from its environment, 
decoherence poses a severe challenge to the application 
of quantum technologies such as quantum information 
and communication [2]. The effect of the environment 
is not necessarily detrimental. For example, it can be 
used for robust implementation of quantum processes [3] . 
Whether one wants to fight or exploit the effect of the 
environment, a profound understanding of the system- 
environment interaction and the ensuing open quantum 
system dynamics is indispensable. 

The challenge for a rigorous treatment of open quan- 
tum system dynamics is due to the system and envi- 
ronmental degrees of freedom (DOF) becoming entan- 
gled because of their interaction. This entanglement is 
neglected when invoking the Markov approximation [T]. 
Non-Markovian effects are abundant, in particular in the 
condensed phase. Examples of current interest include 
solid state devices that interact with defects modelled 
as two-level fluctuators [H [5] or excitons of the light- 
harvesting complex immersed by proteins While un- 
der specific assumptions, theoretical methods to tackle 
non-Markovian environments can be derived [IH91 , a nu- 
merically feasible and generally applicable method to 
study the dynamics and control of open quantum sys- 
tem is lacking. 

A complete microscopic quantum description of both 
system and environment, 



duced dynamics of the system. 
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accounts for any type of system-environment interaction. 
Using Eq. ([T]) is, however, prohibitive due to the expo- 
nential scaling of the Hilbert space dimension with the 
number of DOF. Typically a complete description is also 
not needed. In particular, if one is interested in the re- 
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at finite times, for example in femtosecond photochem- 
istry, the complete description can faithfully be replaced 
by a 'surrogate Hamiltonian' [TUHT2]. This is due to the 
fact that at finite times, the system can only resolve a 
finite part of its environment. 

Here we combine the concepts of the 'surrogate Hamil- 
tonian' and the Lieb-Robinson bound [131 HI] to devise 
a general approach to model decoherence starting from 
a microscopic description of the system interacting with 
its environment. Our existential proof yields an efficient 
algorithm to calculate the reduced dynamics. We use the 
Lieb-Robinson bound to translate the inherent locality of 
quantum dynamics into a quantitative estimate for the 
information propagation speed |15j . The notion of an 
effective light cone can thus be used to set up a dynam- 
ical renormalization procedure for the generator of the 
'surrogate' evolution. 

Graph representation of the full Hamiltonian Assum- 
ing the interactions between system and environment to 
be bilinear and bounded and the environment to be com- 
prised of discrete DOF, the Hamiltonian ([I]) can be de- 
fined on a generic lattice in arbitrary dimensions, 



(3) 



i<j=l 



with Ns system and Nb environmental (or 'bath') DOF 
{Nb — >■ oo). We will later generahze our results to con- 
tinuous environments, e.g., comprised of harmonic oscil- 
lators, and unbounded interactions. In Eq. ([s]), each de- 
gree of freedom is defined on a finite-dimensional Hilbert 
space. Hi, with iV™* (N^*) denoting the number of those 
system (environment) DOF that interact with the envi- 
ronment (system) . The interactions can be expressed 



2 



in terms of operators Oi belonging to the algebra BCHi) 
of bounded linear operators on Hi, 



E 

fj.=0 



dim(B(Hj))-l 



v=0 



where the diagonal terms represent on-site interactions 
and \ Ji^\<oo. The Hamiltonian (Isj) defines a graph, i.e., 
an ordered pair G={T,E) of nodes and edges. The set 
of nodes F — {Ns + Nb} consists of all the system plus 
environment DOF and hence can have an infinite number 
of elements. The set of edges E — {Jij ^ 0} represents 

I 1/2 

. The bare 



all non-zero couplings Jij 
structure of the graph is encoded in the adjacency matrix 
A = A{G) whose entries Aij {Aij — 0, 1) are the number 
of edges connecting i and j. 

Regularity of the graph To apply a renormalization 
procedure to Hamiltonian ([3|, the associated dynam- 
ics has to be quasi-local. This corresponds to the 
graph being regular. Regularity of the graph G is 
equivalent to the existence of a non-increasing func- 
tion F, F{d{i,j)) : [0, oo) [0, oo), acting on a met- 
ric d{i,j) defined on G, that is uniformly integrable, 
= supjgp ^^gp F(c?(i, j)) < oo, and obeys the con- 
volution property, sup, E^gr ■^('^(*' ^))^('^(^' i)) ^ 
GF(d{i,j)) with C<oo. The proof requires us to assume 
that the graph, although infinite, is locally finite. Defin- 
ing the connectivity Ci of the ith node as the number of 
edges emanating from it, local finiteness implies Ci < oo 
for each i. The connectivity encodes both the dimen- 
sionality and the range of interactions on the physical 
lattice which therefore have to be finite. Local finite- 
ness of the graph is necessary and sufficient for the ad- 
jacency matrix A to be bounded in the operator norm 
on G, ||A|| < sup.gp Ci = supjgp J^'jLi ^ij = c US]. The 
boundedness of A ensures boundedness of the coupling 
matrix J. We can thus define the metric on G as the 
shortest path connecting two nodes, 

d{i,j) :=min{neNo : ^ 0} , 

which does not necessarily coincide with the distance be- 
tween sites on the physical lattice. A walk of length 
n from node i to j can then be defined as a sequence 
T^nihj) = [i — *o,ii,''' , «n = j] of n adjacent nodes. 
Defining J = J / || J||, the weight of the walk 7r„(j, j) be- 
comes ■w{iTn{i,j)) = XYkZi Jik,ik+i^ "^ith the weight of 
the zero-length walk set equal to 1. The set n„(i,j) 
of all paths of length n between nodes i and j contains 
|nn(*, j)| = elements and weights w(n„(i,j)) = 

[J"]y. The weight of the shortest path(s) between nodes 
i and j is [J'''*'"'']ij and can be used to define the desired 
non-increasing function F, 
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with r > 1 a damping factor penalizing longer paths. 
It follows from the boundedness of J that Fr{d{i,j)) 
satisfies the uniform integrability condition, \\Fr\\ < 



< 



1_ 
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For a given 
quantifies the maximum influence of the 
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damping r, \\Fr 
graph on a single node and vice versa [25] • Fr{d{i,j)) 
clearly satisfies the convolution property with convolu- 
tion constant C equal to one. 

Lieb-Robinson bound To assess the quasi-locality of 
the dynamics generated by Hamiltonian ([s]), consider a 
finite subset X of G, surrounding the system nodes S, 
together with the truncation of H on this subset, Hx = 
V X^i jex ^^j- Denoting the boundary of the region 
S (containing the system nodes) by dS and that of the 
region X (containing the environment nodes closest to S) 
by dX, the graph distance between dS and dX defines 
the radius of X, i?| = d{dS,dX). The boundary dX 
consists of those elements of X that interact with DOF 
outside of X. Approximating Eq. ([3| by the truncated 
Hamiltonian Hx gives rise, in the Hcisenberg picture, to 
an approximated time evolution, Og{t), for any system 
operator Os ■ The error introduced by the approximation 
is bounded. 



Osit)-0^{t) <2 Os \\Fr\\\dS\e 



(5) 



where Wr=2||l>|jre and ||<i>||r=supjjgp \\<^ij\\Fr{d{i, j))^^ . 
The difference is taken in the operator norm and \dS\ de- 
notes the number of elements of dS. Equation ^ states 
the Lieb-Robinson bound (1^15j for open quantum sys- 
tem dynamics (see the supplementary material for de- 
tails) , reflecting the approximate locality of the dynamics 
generated by H. At any fixed time, the full generator can 
therefore be approximated by a truncated one. The ap- 
proximation improves exponentially with increasing ra- 
dius, Rx, of the effective support, i.e., the contribution 
of environmental DOF outside of the effective light cone 
of S to the system dynamics is exponentially small. The 
bound ([5| is fully general since it relies only on the cut- 
off imposed by the discrete graph structure [SO]- For a 
prespecified error ei, Eq. ([s]) deflnes a time window. 



t e 



:{0,(ln(G, 
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with Cr = ei (2 ||$||^ IIF^II \diS\y^ E (0, 1). The value 
of r can be set to 2 by maximization of t^s and the 
subscript r = 2 will be dropped below. Within the time 
window ([6]), the full generator, H, can be replaced by 
a surrogate one, Hx- The bound i^s accounts for the 
worst-case relation between interaction time and radius 
of the effective support. Therefore, the full environment 
is needed at most in the limit of infinite interaction time. 

Truncation of the environment: dynamical renormal- 
ization group The ability to truncate H for finite times 
provides the basis for a renormalization procedure to de- 
termine, given t and ei, the 'surrogate Hamiltonian' Hx- 
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Reordering the environmental DOF according to their 
graph distance d from S in an onion-like fashion, the full 
Hamiltonian Q is rewritten, 



H 



d=Q 



{hd + hd^d+i) 



(7) 



where hd accounts for the interactions between all DOF 
at a distance d from 5, i.e., within the layer dXd — 
Xd \ Xd^i, while interactions between two successive 
layers, dXd and dXd+i are grouped within hd.d+i- The 
lowest order terms include the system DOF, ho — Hs, 
fiQ^i = HsB- After the reordering, the DOF in Xi are 
labelled by indices i — (1, • • • ,A^^"*), those in 8X2 by 
i = (iV™*-|-l, . . . , Ni;^*+\dX2\), and those in the dth layer 

aXdbyz = E^t!|5^jl + l,--- ,E,til^^jl)- Truncat- 
ing the sum in Eq. ([7| yields a set 




H 



k=l 



fe) ) 



(8) 



which approaches H in the limit d— > 00. Correspond- 
ingly, the bound Q defines a sequence of increasing times 
{td}T=o with td = max [0, (ln(C) + d) /v]. 

Then the following dynamical renormalization group 
transformation for the generator holds 

t — td =^ Hxd-i^Hxa= Hxd-i+hd + hd-i^d 



VdeNo, i e [td,td+i), (9) 

which ensures that the error in the simulation is kept 
fixed by retaining, within each time interval, the corre- 
sponding minimal set of environmental DOF, cf. Fig. [T] 
Clearly, the number of environmental DOF required at 
each step is model-dependent. 

Efficient simulation For a given t, the computational 
complexity of simulating the open quantum system dy- 
namics can be further reduced by applying the Suzuki- 
Trotter expansion: For td < t < td+i, the effective prop- 
agator defined by Hx^ can be approximated as 




(10) 



where At — t/nid and the product runs over all the 
interacting pairs {i,j} belonging to Xd taken in arbi- 
trary order. The generator Hx^ contains Kd < \Xd\'^ 
terms. The error introduced by approximating e"'^-^'*'^* 
within each At by a product of Kd terms is of the 
order £2 < \c%Kl{MY with c| = sup{,j}g^^_^ 

[T71[THj. A prespecified error 62/2 for the whole time t is 
achieved by taking = c\t^K^/e2 Trotter steps, i.e.. 



FIG. 1: (color online) Pictorial representation of the renor- 
malization transformation (Eq. (|9|): When t reaches td, the 
effective generator Hx^-i is updated by adding the interac- 
tions with the environmental DOF in the dth layer and the 
new local terms. 



At = €2/{tc\K'^). The characteristics of the microscopic 
description (|3|, dimensionality and interaction range, de- 
termine only Kd- An increase in either one of them there- 
fore increases only the number of Trotter steps md but 
not the size of the operators. For any time t the full dy- 
namics of system and environment, \^Xa{t)), can thus be 
determined according to Eqs. ([9| and (|lO|. The state of 
the system, ps{t), is obtained as tr^ [|^'Xd(i))(*Xd(OI] 
and differs from the true one evolved according to Eq. ^ 
by at most ei + £2/2. 

Extension to continuous environments The above 
renormalization scheme generalizes to both fully con- 
nected and unbounded interactions. Indeed two stan- 
dard models of decoherence, the central spin model and 
the harmonic bath ill, fall into these categories. As a 
benchmark example, consider a system interacting with 
a continuous bosonic environment, 



H^Hs + Os / h{e){de + dl)de + / ^(e)^^^^*, 
Jo Jo 

(11) 

where d^ (a| ) denotes the annihilator (creator) for a bath 
mode of energy e, and Wc is a cutoff frequency. The bath 
dispersion g{e) and the system-bath coupling h{e) are re- 
lated to the spectral density J(e) = '^''^^^-d^h^ig^^i^)) 
[in] . By a proper choice of basis using orthogonal poly- 



nomials, every Hamiltonian of the form (11) is unitarily 
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equivalent to a semi-infinite chain [501 HI] , 

CO 

H = Hs + Jods{ba + bl) + ^uJny^J)n+tn{blK+i+h.c.). 

(12) 

This establishes local finiteness for continuous environ- 
ments and holds also for discrete Hamiltonians of the 



functional form (11). Since (11) is quadratic in the 



bath operators, a Lieb-Robinson bound holds even for 
unbounded interactions |22| . In the supplementary ma- 
terial, we present a direct derivation of the bound. It is 
based on Hb being quadratic such that the free evolution 
of the environmental DOF is mapped into that of wave 
functions / S P{T) on the graph, = Ylt^o fnk{t)ik, 
with fnk{t) — [e~*'^^*]„fe and Jb the bath coupling ma- 
trix. The wave functions are bounded, 

<^;(|n-fc|)e-(l"-'=l-"'-*), (13) 

mm{JB.tj)/ \\Jb\\- 



with = er II Jb|| /J^„^ b^ Jmin,B 
p: 

generator becomes 



In the interaction picture, 6' (t) = e^HBtQ^-iHst ^ 



H' (t) = Hs+HsB (t) - Hs+JoOs V (fok it)h+fSk (i)^l 



fe=0 

(14) 

and, using the Suzuki- Trotter expansion, the cor- 
responding propagator is written as V^{t) ~ 
lira„,^ooUT=o^(tvh + ^t)^ with y(tj,tj + At) = 
^-iHsAt^-iHsBijAt) ^ As shown in the supplementary 



material, the bound (13) yields the following estimate 



for the accuracy of the effective propagator defined 
by H^{t) accounting only for the first d bosons in the 
chain dltl), 



oo „t 



k=d 

< \\F;'\\{2vr) -e 



2||(2t;,)-ie-2('^--^*). (15) 

This is a Lieb-Robinson bound analogous to Eq. ^ 
with r = 2 again and serves as the starting point to set 
up the renormalization transformation ([9| for the har- 
monic oscillator bath. Since the interactions are nearest- 
neighbour in a one dimensional chain, no reordering of 
the environmental DOF is required prior to renormaliza- 
tion. 

Conclusions We have shown that the reduced dynam- 
ics of a generic open quantum system can be determined 
efficiently and with arbitrary precision. The only as- 
sumption to be fulfilled is that each degree of freedom 
interacts with a finite number of other DOF. This holds 
for interactions that decay with distance or, in particu- 
lar for continuous systems, if a finite energy cut-off ex- 
ists. The computational effort scales at most polynomi- 
ally in time and size of the effective support of the en- 
vironment. The memory required to store the surrogate 



Hilbert space vector l^'x^) scales exponentially with the 
number of DOF in the effective light cone. This scaling 
can be made polynomial by restriction to excitation sub- 
spaces [TT]. Our result can be extended to account for 
arbitrary time-dependent bounded interactions applying 
available generalizations of the Suzuki- Trotter decompo- 
sition While generalizations of the Lieb-Robinson 
bound to dissipative dynamics exist [24l |25] , they have 
never been applied to the full system-bath evolution. Our 
work thus provides a novel and very natural application 
of the bound. 

The bipartition of system and environment intrinsic to 
our renormalization scheme makes it more powerful for 
modelhng decoherence than t-DMRG This is 

due to the existence of an effective light cone that can 
be fully exploited in terms of an optimized truncation 
of the environmental DOF if one only cares about the 
dynamics of a small subsystem. Our renormalization 
scheme allows to go beyond one-dimensional lattices with 
nearest-neighbor interactions while keeping the quality 
of the approximation throughout time. By resorting to 
hyper-graph theory, our proof is extended to arbitrary 
fc-local interactions (as shown in the supplementary ma- 
terial). Our method thus opens the way to study a huge 
range of open quantum systems, including those in non- 
Markovian environments, and their control. 
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action Q0C4QIP Grant No. PIEF-GA-2009-254174 is 
gratefully acknowledged. 



[1] H.-P-Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford Univ. Press, 2007). 

[2] M. A. Nielsen and I. L. Chuang, Quantum Computa- 
tion and Quantum Information (Cambridge Univ. Press, 
2000). 

[3] J. T. Barreiro, M. Miiller, P. Schindler, D. Nigg, T. Monz, 

M. Chwalla, M. Hennrich, C. F. Roos, P. Zoller, and 

R. Blatt, Nature 470, 486 (2011). 
[4] G. de Lange, Z. H. Wang, D. Riste, V. V. Dobrovitski, 

and R. Hanson, Science 330, 60 (2010). 
[5] H. Bluhm, S. Foletti, I. Neder, M. Rudner, D. Mahalu, 

V. Umansky, and A. Yacoby, Nature Phys. 7, 109 (2010). 
[6] M. Sarovar, A. Ishizaki, G. R. Fleming, and K. B. Wha- 

ley, Nature Phys. 6, 462 (2010). 
[7] J. Piilo, S. Maniscalco, K. Harkonen, and K.-A. Suomi- 

nen, Phys. Rev. Lett. 100, 180402 (2008). 
[8] H.-P. Breuer and B. Vacchini, Phys. Rev. Lett. 101, 

140402 (2008). 

[9] W. Koch, F. Grofimann, J. T. Stockburger, and J. Anker- 
hold, Phys. Rev. Lett. 100, 230402 (2008). 
[10] R. Baer and R. Kosloff, J. Chem. Phys. 106, 8862 (1997). 
[11] C. P. Koch, T. Kliiner, and R. Kosloff, J. Chem. Phys. 

116, 7983 (2002). 
[12] C. P. Koch, T. Kliiner, H.-J. Freund, and R. Kosloff, 



5 

Phys. Rev. Lett. 90, 117601 (2003). 
[13] E. H. Lieb and D. W. Robinson, Commun. Math. Phys. 
28, 251 (1972). 

[14] B. Nachtergaele and R. Sims, in New Trends in Math- 
ematical Physics. Selected contributions of the XVth In- 
ternational Congress on Mathematical Physics, edited by 
V. Sidoravicius (Springer, 2009), pp. 591-614. 

[15] S. Bravyi, M. B. Hastings, and F. Verstraete, Phys. Rev. 
Lett. 97, 050401 (2006). 

[16] B. Mohar and W. Woess, Bull. London Math. Soc. 21, 
209 (1989). 

[17] H. F. Trotter, Proc. Am. Math. Soc. 10, 545 (1959). 

[18] M. Suzuki, Proc. Japan. Acad. 69, 161 (1993). 
[19] R. Bulla, T. A. Costi, and T. Pruschke, Rev. Mod. Phys. 
80, 395 (2008). 

[20] J. Prior, A. W. Chin, S. F. Huelga, and M. B. Plenio, 

Phys. Rev. Lett. 105, 050404 (2010). 
[21] A. W. Chin, Angel Rivas, S. F. Huelga, and M. B. Plenio, 

J. Math. Phys. 51, 092109 (2010). 
[22] B. Nachtergaele, B. Schlein, R. Sims, S. Starr, and V. Za- 

grebnov. Rev. Math. Phys. 22, 207 (2010). 
[23] D. Poulin, A. Qarry, R. Somma, and F. Verstraete, Phys. 

Rev. Lett^ 106, 170501 (2011). 
[24] D. Poulin, Phys. Rev. Lett. 104, 190401 (2010). 
[25] M. B. Hastings, Phys. Rev. Lett. 93, 140402 (2004). 
[26] G. Vidal, Phys. Rev. Lett. 93, 040502 (2004). 
[27] S. R. White and A. E. Feiguin, Phys. Rev. Lett. 93, 

076401 (2004). 

[28] G. De Chiara, M. Rizzi, D. Rossini, and S. Montangero, 

J. Comp. Theor. Nanosci. 5, 1277 (2008). 
[29] In graph theory, \\Fr\\ represents the maximum Katz cen- 

trahty of the nodes [L. Katz, Psychometrika. 18, 39-43 

(1953)]. 

[30] In some cases tighter model-dependent bounds can be de- 
rived [Burrell and Osborne, Phys. Rev. Lett. 99, 167201 
(2007)] 



6 



Supplementary Material for 
Approximating open quantum system dynamics in a controlled and efficient way: 
A microscopic approach to decoherence 

Giulia Gualdi and Christiane P. Koch 
Theoretische Physik, Universitdt Kassel, Heinrich-Plett-Str. ^0, D-34132 Kassel, Germany 

I 



QUASI-LOCALITY OF QUANTUM DYNAMICS 

In this section, the proof of Eq. (5) in the main text is 
sketched for the sake of self-containedness, following Refs. 
[m [TS] . Two steps are required: (i) proof of the Lieb- 
Robinson bound in its original formulation as a bound 
on the commutator of two observables taken at different 
times |14| , and (ii) demonstration that this bound implies 
that the time-evolution of any observable can be well 
approximated by an operator acting only on degrees of 
freedom within its light cone [T5] . 



Lieb-Robinson bound 

Consider a set F of degrees of freedom, labeled by 
the index z, each corresponding to a finite Hilbert space 
T-Li. The full Hilbert space of the system is then T-Ly — 
0jgP 'Hi- For each degree of freedom, the observables be- 
long to the algebra B{'Hi) of bounded operators on Hi- 
Given a finite subset X C T, the algebra of observables 
over X is defined as Ax = ®i(zx^^i)- By identifying 
an operator A G Ax with A(E)Ir\x G ■^r, it then follows 
that Ax C At- If the set F is infinite, Ar is defined by 
the inductive Hmit = ^xerAx- The support of an 
observable A G Ar is defined as the minimal set X C T 
for which A can be written A = A' (E) ir\x with A' E Ax ■ 
A local observable is therefore an observable whose sup- 
port is a finite subset of F. 

The interactions among elements of F can in general be 
defined as a map from the set of subsets of F, {X}xcr, 
to the algebra Ar, $ : {X}xcr Ar, which obeys the 
following properties 

. ^X)eAx; 

A generic Hamiltonian on F can then be written 



H 



(16) 



and, for any observable A € Ar, the dynamics generated 
by H is expressed as 



rt{A) 



(17) 



Intuitively, the spread of the support of an observable 
A due to its dynamical evolution should depend on the 
surface area rather than the volume of its support. The 
surface of a set A C F is the subset of F having non-zero 
intersection both with X and its complement A^, 

S{X) = {Z CT : ZnX ^0 and Z n A V 0} . 

The boundary of a set dX is then defined as the subset 
of X whose elements interact with degrees of freedom 
outside of A, 

dX ^ {x e X : 3Z e S{X) with x G Z and $(^) ^ 0} . 

Finally, given a metric d on F, the distance between 
two sets A C F and r C F is defined by d{X,Y) = 
minxex.yeY dix,y). Equipped with these definitions the 
following holds [13] : 

Lieb-Robinson bound. Given a set F and a met- 
ric d on it, if there exists a non-increasing function 
F : [0, oo) — > [0, oo) with the following properties 

1. uniform integrability: 



2. convolution: 
3Ce 



-<oo: supV»^M^<c, 



then for any fi G ]R\j{0} and any pair of local observables 
Ae Ax and B e Ay with d{X, Y) > 0, 



xer 



rt{A),B 



< 



mum 

Cu 



|F||min{|9A|,|9r|}e" 



-t,{d(x,Y)-2\m,,c^\t\) 



(18) 



with 
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1. m\ 



2. the convolution constant corresponding to F^{d{i^ j)) — e '^'^*^*'^^_F(c?(i, j)). 



Equation ( 18 ) states that the support of a given ob- The first term on the right hand side of Eq. ( 19 ) is norm 



servable spreads at a velocity v with 



V < inf 

Ai>0 



To prove Eq. ( 18 ) the foUowing quantity needs to be es- 
timated 



CB{X:t) = sup 
AeAx 



rtiA),B 



< 



assummg that T,tedxT,jeY ^W'^^j)) 
X^jear Siex ^(^(*'-?)) (^^ ^^^^ assumption is not 
fulfilled the proof below needs to be carried out for 



T-t{B),A 



It is useful to isolate in Hamiltonian ( 16 ) 



loc I ^<E> 



the on-site terms hi, 
The dynamics generated by 



zex 



leaves the algebra Ax of observables on X invariant, i.e., 
the propagated operator defined for A E Ar, 



loc 



wiU again be in Ax, tI°''{A) e Ax, if A e Ax- Using 
the property of the commutator, one can then write 

fit)^[n{A),B]^[nirLriA)),B] . 

Differentiating f{t) with respect to time yields the fol- 
lowing differential equation. 



fit) - ^ E [rtmZ)J{t) 



zes{x) 
'i 



E [n(riT(i)), .(19) 

zes{x) 



preserving. Using the properties of first order inhomoge- 
nous differential equations, one obtains 



fm< [A,B] +2\\A\\ E f [rsmZ)),B] 



zes(x) ■ 



ds . 



(20) 

Due to inequality (20) and the fact that t/°'^(-) is norm 
preserving, 



zes{x) 



CB{X,t)<CBiX,0) + 2 V U{Z) / CB{Z,s)ds, 

(21) 

where Cb{X; 0) = since we assumed d{X, Y) > 0. Iter- 
ating the same argument used in (21) for the integrated 
function yields 



CB(X,t)<2||i?||£M^ 



n=0 



with coefficients a„ (n > 1) 



E E 



E Sy{Z^)l[U{Zi) 



ZiGS(X) Z2eS(Zi) Z„GS(Z„_i) 



1=1 



(22) 

where 5y{-) = 1 if F n • 7^ and zero otherwise. To 
put an upper bound on the coefhcients a„ requires exis- 
tence of the function F{d{i,j)). After some algebra and 
employing the convolution property of F, the following 
estimate is obtained, 



an<mrc-' E EW'-?'))' (23) 

where |||.|| = sup^^r ll^-WIKW, j))-'- In- 

serting (23) into (21), one obtains 



[n(i),S]||<M^(e2C||*l||,_i) y: EW,J)) 



(24) 



which is the first form of the Lieb-Robinson bound. Note that if the supports of the observables A and B have a 
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non-empty intersection, i.e. d{X, Y) = 0, then one only 
needs to substitute (e^^H*!! 1*1 -1) by e^^ll*!! 1*1 in 1241. To 



recast (24) in the form of (18 1, it suffices to note that the 



properties of F{d{i,j)) are obeyed also by any function 
of the form F^{d{i,j)) = e-'"^(*'^)i^(d(z, j)) with > 0. 



Indeed, ( 18 ) is obtained by replacing F by in Eq. ( 22 ) 
and, correspondingly, $ by 



I'll, = sup mx)\\{FM^-j))r' 



and noting that 



<e-'^'^(^'^)||i^||min{|a*X|,|a*r|} . 



Approximate locality of quantum dynamics 



The second step consists in showing that the 



bound (18) implies quasi-locality of the dynamics gen- 



erated by H of Eq. (16). In other words, the dynamical 



evolution Tt{A) of a local operator A with support in 
X C r can be well approximated by an operator acting 
only on the degrees of freedom within the effective light 
cone of A. As shown in Ref. [TT, this operator can then 
be identified with the time evolution of A generated by 



those terms in (161 that are supported on the effective 



light cone. In the following section, we will establish an 
analogous result for continuous systems. 

Denoting by C F the set of degrees of freedom hav- 
ing at least distance d from X, the (normalized) partial 
trace of Tf [A) over Y'^ is an operator acting on 'Hy 3 Tlx 
defined by 



(n(i)>i 



tryc(Tt(i))0lyc 
trye(Iyo) 



e Ay 



(25) 



It is convenient to represent the partial trace as an invari- 
ant integration over the group of unitary transformations 
Uiy^) acting on the Hilbert space corresponding to Y"^, 

(rt(i))y= / {lY®U^)A{lY®U)d^l{U) (26) 
with normalized Haar measure dfi{U), 



d^{U) — di'm{'HY<') ■ 



U{Y-=) 



Equations ( 25 ) and ( 26 1 define indeed the same opera- 



tor since both expressions, starting from Tt{A), uniquely 
define an operator which commutes with all operators of 



the form ly (g) O. Using (181 and (26) one obtains 



< 



< 



u{y= 

c,. 



Tt{A),U] dfi{U) 



which is precisely a Lieb-Robinson bound of the form of 
Eq. (4) of the main text. 



LIEB-ROBINSON BOUND FOR HARMONIC 
ENVIRONMENTS 

This section is devoted to the proof of Eq. (13) in 
the main text. In Ref. [T] a Lieb-Robinson bound for 
quadratic Hamiltonians with unbounded interactions has 
been obtained. Here we present an alternative simpler 
proof of the bound tailored to the kind of quadratic 
Hamiltonians that is used to model decoherence due to 
an environment consisting of harmonic oscillators. 

The goal is to obtain an explicit estimate of the 
quasi-locality of the dynamics generated by the follow- 
ing Hamiltonian 



Hs + HsB +Hb^Hs + JaOs {ba + b, 



OO 

■J2^nbibn+t„ (blX+i+h-c.'^ ,(27) 



n=0 



where Hb can be written as Hb = bJsb with b = 
(6o, • • • ,bn) and Jb the coupling matrix of the environ- 
mental degrees of freedom with entries a;„ and Then 
the normal mode representation of Hb, 



ai , 

can simply be obtained by diagonalizing Jb, 

Xi = [UJbU^)i ,k^Y. ' = E ^Ikk 

I k 

Since the normal mode free evolution is given by 



the operators in the original basis (27 1 evolve according 
to 



n 

= ^[C/diag(e-'^'*)C/t]^J„ 

n 

= ^(7A„fe(<)6„-f Cfe(i)St), (28) 
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where A/ denote the eigenvalues of Jb- Hence, for any 
initial linear superposition of bosonic operators, 

Kf) = Y.^fkbk + rA) , 

k 

with / = (/i, ■ ■ • , /oo) e /^(r), one can write 
TMf)) = KTtf) = 6 ([C/diag(e-'^'*)C/t] /) 

k,n 

= J2 fr^kit)bn + fnk{t)bl (29) 



In the interaction picture defined by 

d'{t)=e'^^'6e-'"''\ (33) 
the generator of the evolution of the states becomes 



H'it) - Hs + JodsJ2(fokit)bk + fokit)bl) 



k=0 



Hs + HsBit). 



(34) 



Employing the Suzuki-Trotter expansion [I7l[TH], the cor- 
responding propagator can be written 



It follows from Eq. (29) that the existence of normal 



modes allows for mapping the dynamical evolution of op- 
erators on the Hilbert space Jir into that of wave func- 
tions on r. 

The task is now to bound each \fnk{t)\- To do so we 
set 



B hi 



j.d{i,j) ' 

where Jb = r > 1 and 

d(i,j) = min{neNo:[A"],j^O} 

with A the adjacency matrix of Jb and \\Jb\\ < ||^|| < 2 
since A is tridiagonal. Defining 



Fr{m)= sup Fr{d{i,j), 

i,j-d{i,j)=m 



we can estimate 



\ min J 



(30) 



(31) 



with Jmm = mm(^ij);jg ^.^oJB,ij- Using Eqs. (|30j) and 
(31 1, the following bound is obtained 



\fnkit)\ 



-iJst] 



.fe/fc(0)| 



■[Jmnkfkio) 



< 



E 

E —^ \[JB]nke''''^Fr{m)e^'^{Fr{m))- 
^ — ' to! 



< e-H"-felF,(|n-A;|) V — (||Jb|| </.,)" 

^ — ^ m! 

in—\n—k\ 

< F^(|n-/c|)e^''(l"-'=l-<*) (32) 



with ^ > and u^' = \\Jb\\ |e'^0r/M- Equation (32) is 
analogous to ( 18 ) and provides the key to estimate the 



approximate locality of the dynamics. The freedom in 
the choice of fi can be removed by minimizing which 
yields = 1. 



V'{t) = lim ]T e-''^^^*e-^^^^ 

m. — ^na J- J- 



O'At) 



3=0 



lim TT V{t^,t^ + At) 



(35) 



where At — t/m. The propagator associated to the trun- 
cated generator, 



HiB{t) = Jods J2 {fokit)h + fSk{t)bl) , (36) 

k=0 

is expressed analogously. Setting 

oo 

HiBit) = JoOs {fok{t)bk + fSkitPk) , 

k=d+l 

the error made by employing the truncated generator 
within an infinitesimal time interval can be estimated 



as 



\\V{t„tj+At)~Vdit,,t,+At)\\ 

= j|p-iHsB(jAt) _ g-iffsB(iAt)g-i_f/|B(jAt)| 



(37) 



Recognizing that e^^sBU^*-) and e ^fissUAt) ^^^^ 
Weyl operators since they are of the form 

W{f)^exp[ib{f) + b\f)] 

with / G ^■^(r) and hence obey the Weyl relation 

W{f)Wig)=W{f + g)e-''"'^<f-<^^^ 

with (•, •) the inner product on ^^(r), one obtains the 
following estimate for Eq. ( 37 ) 



[e»-f^sB(iAt)^g-iffsBOAt)j 



<|(/o(At),5o(At))| (38) 
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with /o(At) = {/ofc(Ai)}r=o' 5o(Ai) = {/ofc(AO}r=d+i- 
Together with Eq. ([32]), Eq. ([38| yields the following es- 
timate 

< ||F2||e-2M(fe-^.At) ^ 

Taking the limit m — >■ oo, one obtains Eq. (15) of the 
main text, 



the node, c; 



If all nodes have the same con- 



nectivity c, then the hypergraph is c-regular, in analogy 
to the notion of homogenous graphs. Analogously to the 
adjacency matrix of a graph. A''' is bounded if and only 
if the hypergraph is locally finite, i.e., Ci < oo for each 
i e r, and we obtain \\A^\\ < c = sup.gr Ejli "^ij- The 
metric on the hypergraph can then be defined in terms 
of the shortest path connecting two nodes, 

d{t,j) :=min{neNo : [(A'')"],^ ^ 0} . 



\\V\t)^Viit)\\< / e(.)ds<||F,2||(2«.)-ie-2(^-'"-*) 



GENERALIZATION TO fc-LOCAL 
INTERACTIONS 

We now show how the renormalization procedure gen- 
eralizes to Hamiltonians exhibiting arbitrary fc-local in- 
teractions. Since the Lieb-Robinson bound is fully gen- 
eral, the only step required is to define an appropriate 
metric d on the set of degrees of freedom F and a non- 
increasing function F{d{i,j)) obeying both uniform inte- 
grability and the convolution property. 

Consider a generic /c-local Hamiltonian 

(39) 

where, without loss of generality, we have assumed all 
to be A:-local (instead of at most fc-local). 



E 

{mi},--- ,{a's 



} 



3 = i 



(40) 



In Eq. (40), {/Zi} denotes that each index /i^ varies be- 



tween and the dimension of the finite Hilbert space of 
the {ji)th degree of freedom, dim{'Hj.), and Oj- stands for 
the tensor product of the operator Oj. e BCHj- ) and the 
identity over the other degrees of freedom. The Hamil- 
tonian (39) defines a hypergraph, i.e., an ordered pair 
Q = {T,E) where the set of nodes T represents the de- 
grees of freedom and E is the set of non-empty subsets 
of r, called hyperedges or links, for which Jf^'...' 7^ 0. 
Since we have assumed all the interactions to be fc- 
local, all hyperedges have size fc and the hypergraph is 
fc-uniform. A graph can therefore be regarded as a 2- 
uniform hypergraph. The adjacency matrix A^^ of an 
hypergraph Q is defined as the matrix whose entries A^,^ 
correspond to the number of hyperedges containing both 
degrees of freedom i and j [2 . The connectivity of a 
node Ci is given by the number of hyperedges involving 



One choice for F{d{i,j)) is then 



Fr{dii,j)) 



j.d{i,j) 



with A'^ = With these definitions we can 

proceed analogously to the case of bilinear interactions 
and define the dynamical renormalization group trans- 
formation (Eq. (8) in the main text). Once an effec- 
tive propagator for td < t < td+i has been obtained, it 
can be approximated using the Suzuki- Trotter decompo- 
sition [TTlfTSI. 



n 



where the product is over all time intervals A^ — t/md 
and all fc-tuples belonging to Xd. We find that for a 
prespecified error £2/2 the effective propagator for time 
t can be approximated as a product of Ud — (aptKd)^ /e2 



terms where ci 



and Kd < \Xd\ the num- 



ber of fc-local unitaries. One can then use the Solovay- 
Kitaev algorithm [3l [23j to further decompose each fc- 
unitary into a product of one- and two-body unitaries 
chosen from a suitable set. To achieve an accuracy e for 
each fc-unitary transformation, nsK = ftlog2(e~^) oper- 
ations are required with a, b constants. If one chooses 
e = £2/ (2nd), the effective propagator is simulated with 
an accuracy £2 employing — and\og2{nd/e2) one- and 
two-body unitaries, i.e., with a computational effort that 
scales polynomially in time and size \Xd\ of the effective 
light cone. 
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